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Abstract 

We study an effect of the off-shellness of the quark and antiquark inside a heavy 
quarkonium system on IR renormalons contained in the perturbative computations of the 
quarkonium energy levels. We demonstrate that, when the off-shellness pq — niQ ~ a'^mq 
is larger than mg Aqcd, renormalons in the energy levels as calculated in perturbative 
QCD are suppressed by a factor Aqcd / ("I'tt-q) as compared to those in 2mpoie + Vqcd (5") • 
In this case the residual C(Aqqj-,) renormalon has the same dimension as that of the 
leading gluon-condensate contribution. 
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Recently, there has been renewed interest in infrared (IR) renormalonsQ in the context of 
perturbative computations of the physical quantities of the heavy quarkonium states, within 
the framework of effective theories such as Non-Relativistic QCD (NRQCD) P| and potential- 
NRQCD (pNRQCD) 0]. This is because, proper identification of renormalons and realization 
of their cancellations in these perturbative computations have led to significant improvements 
in the determinations of the heavy quark masses m^, rrib and nic 0, ||, ||, ^ ||, 0, as well 
as in the perturbative computations of the heavy quarkonium spectra 0, ^ and the static 
QCD potential [1^, |11|. We have learned that through a deeper understanding of renormalons 
in these systems we may achieve more accurate perturbative predictions, and consequently 
we may gain more accurate physical picture of the heavy quarkonium states. So far, our 
understanding of IR renormalons in the perturbative computations has been as follows. (See 
also the discussion at the end of the paper.) For simplicity we consider the case where the quark 
Q and antiquark Q have equal masses in this paper. When the total energy of a static QQ pair 
Eiotir) = 2mpoie + Vqcd('^) is expressed in terms of a short- distance mass of Q{Q), where mpoie 
is the pole mass of Q{Q) and Vqcd('^) is the static QCD potential, the leading renormalon 



of order Aqcd contained in 2mpoic [1T2|, |T3[ (the self-energies of Q and Q) is cancelled |jT^, [15 
against that contained in Vqcd('") |16[ (the potential energy between Q and Q). Remaining 
sub-leading renormalons in Etot{r) are of order Aqcd x (^^Aqcd)^ [T^l' "where r is the size 
of the boundstate. Hence, they are smaller than the leading renormalon contributions for a 
heavy quarkonium system, for which tAqcd ^1. As a result, the perturbative expansion of 
the total energy Etot becomes much more convergent when the expansion is expressed in terms 
of a short- distance mass than when expressed in terms of the pole mass. 

In perturbative QCD calculations, the order Aqcd leading renormalons may be identified as 
follows |]I5[. The pole mass and the QCD potential constitute the total energy of a QQ system 



in which Q and Q are on-shell. Dominant contributions from IR regions to these quantities 
can be written (independently of the quark short- distance mass) as 

/■ d^k asi\k\) 

2mp.e~-I/QCD-y^^. (1) 

The integral is dominated by contributions from ~ Aqcd where q;5(|A;|) becomes large. 

The above argument, however, needs to be modified in the case where Q and Q are 
persistently off-shell. In a realistic quarkonium system, Q and Q are off-shell due to the 
binding energy. The off-shellness in a sufficiently heavy quarkonium system is given by 
Pq — ~ 0{alm?). Throughout this paper, m represents the short-distance mass of Q 
and Q, which differs from the pole mass (formally) by an amount of order aim or smaller 
within perturbative QCD, e.g. the Potential-Subtracted mass with an appropriate choice of 
the subtraction scale |jl5| . 



Let us consider the time evolution of a quark-antiquark pair bound in a heavy quarkonium 
state. The potential Vqcd('") generated by a static color source contains a renormalon of 
C(Aqcd)- We may infer that this feature is related to the fact that vacuum polarization by 
gluons is induced around the bare color source, and that an energy of O(Aqcd) is accumulated 
* See [0 for a review of renormalons. 
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Figure 1: A diagram showing typical propagation of Q and Q in a heavy quarkonium state. 



due to the antiscreening effects. When the quark (color source) is static, the antiscreening would 
be completed in a time scale At ~ Aq^d needed for the gluon clouds to surround the quark. 
This picture is consistent with the observation |l5l that the leading renormalons contained 
in the pole mass and the static QCD potential, which are associated with on-shell quarks, 
stem from contributions of instantaneous gluons (in the quark rest frame, and independently 
of gauge choice). Since the on-shell quark is static (stable for infinite time), we may regard the 
time scale At ~ Aq^d as instantaneous comparatively. On the other hand, in a realistic heavy 
quarkonium state, Q and Q are off-shell. This means that Q and Q are not static color sources 
but they keep rescattering with each other at a time interval At ~ l/(a|m). Therefore, if 
the off-shellness is large l/a|m -C 1/Aqcd, Q or Q at rest (in one frame) would be scattered 
before the antiscreening is completed, and the energy accumulated around it would be smaller 
than C»(Aqcd). 

Fig. 1^ shows typical propagation of Q and Q in a heavy quarkonium state. In the leading 
kinematical configuration contributing to the formation of the non-relativistic boundstate, 
gluons exchanged between Q and Q have momenta ~ (9(a|m), \k\ ~ 0{asm). In an 
approximation valid for this configuration, we may replace the gluon vacuum polarization 
as n(A;^) ~ n(— |A;p). Then we may easily integrate over the time-component of the loop 
momenta, and the integral is dominated by contributions from kinematical regions where Q 
and Q are nearly on-shell. On the other hand, it is not obvious whether the same approximation 
is valid also in estimating IR renormalons, since a different kinematical region, — fc^ ~ Aq^j-,, 
contributes to them. In fact, two scales ~ Aqcd and pq — m — pQ/2m ~ a|m compete 
in the quark propagator when Pq — m? ~ a\m? and /c*^, ~ Aqcd, since {pq + k^ — m? ~ 
2m p^Q + k^ — m — [pg + /c)^/2m . Furthermore, when pq — rn^ ~ a|m^ and {pg + k)"^ — 

rn? ~ 0, it follows that k^ ~ a|m, \k\ ~ asm, i.e. k^ ~ a'^m'^, hence in this (nearly on- 
shell) configuration Il{k^) is small and does not contribute to renormalons. In particular, in 
estimating renormalon contributions by power counting in the loop integral of the diagram of 
Fig. D, we should count 

dk° m Aqcd ^2) 



27T {pq + ky — m^ a;|m ' 

if a|m ^ Aqcd- This order estimate is suppressed as compared to the estimate in the on-shell 
approximation of 0{1). According to the power counting (Q), the leading renormalon contained 
in the potential energy between Q and Q is estimated to be of order Aqcd x (AQCD/(a|^)) 
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instead of order Aqcd- Similarly, the leading renormalon contained in the self-energy of Q or 
Q is of order Aqcd x (AQCD/(tt|"^)) and is suppressed in comparison to the renormalon of 
order Aqcd in the pole mass, if we take into account the off-shellness. We will show the power 
countings more explicitly in our discussion below. 

The above argument suggests that, in order to estimate renormalon contributions to the 
heavy quarkonium state accurately in the case a|m 3> Aqcd, we should make an approximation 
valid both in the leading kinematical configuration of the boundstate and in the renormalon 
configuration. Estimates of renormalons within an effective theory may in general differ from 
the renormalons of perturbative QCD, if the effective theory fails to incorporate the off-shellness 
appropriately. Therefore, in this paper, we start from the Bethe-Salpeter (BS) formalism (see 
e.g. ||18|), where one may incorporate (in principle) all the contributions of perturbative QCD, 
in estimating renormalons in the energy levels of boundstate. According to the BS formalism, 
the potential energy Ep^t and the self-energy contributions i?sE in the total energy of the 
boundstate may be expressed, in the rest frame of the boundstate, as 

E^o. = ^{x-K-x) , (3) 



^SE = M - ^ (x ■ [(S^,q)-^(5^,o)"'] ■ x) . (4) 



Here, x?X and K denote the BS wave functions [|I9| and the BS-kernel, respectively; M is 
the boundstate mass (total energy); Sp represents the full propagator of Q or Q. We set the 
momentum of the center of gravity as (M, 0). The dot (■) represents contraction of spinor 
indices and an integral over the relative momentum between Q and Q. Diagrammatically, i^pot 
(£Jse) represents the contributions from the diagrams where the Q and Q lines are connected 
(disconnected), and E^^t + -^^se = M. 

As an example, we estimate renormalons contained in the potential energy iJpot of the 
IS state. Following the standard procedure, we replace the BS-kernel K by the Coulomb 
propagator (in Coulomb gauge) in the large-/3o approximation 



K ^ K^,(k) = -i ^^^^ (7° ® 7») X t (^)" {'og (^) }" ■ (5) 

Here, as{fi) is the coupling constant in the MS scheme, and n is the renormalization scale. 
/3o = 11 — denotes the one-loop coefficient of the beta function with ni flavors; jl = e^^^fi; 
Cp = 4/3 is the color factor. Similarly, it would be desirable to replace the BS wave functions 
x(p),x(p) by those in the large- /5o approximations. Since, however, they cannot be obtained 
simply, in our calculation we substitute the leading-order wave functions in 1/c expansion: 

Xip) - Xfkp) = ( m^ ' P^^- + M ' P^^ ] ^^cAp) X At. 

(6) 
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Here, (j)c,is{p) and M denote the IS* Coulomb wave function and its energy level, respectively: 



2tt {CpdsilA'^)'^ 



2 ' 



p 2 + {CFas{^)m/2f 
A is a spinor matrix representing the spin of the boundstate. 



M = 2m - 



m. 



X 



j=i 



1+7° ^ 1 
V2~ 
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A 



j=o 



1 + 7° 75 1-7° 



(7) 



(8) 



2 V2 2 ' 2 72 2 

with polarization vector e for the vector IS* state. Power counting shows that, as far as the 
dependence on the wave functions x{p)jX{p) is concerned, the order estimates of renormalons 
are determined by the analyticity of by the support in momentum-space of the 

part corresponding to (pcisip)- Therefore, we conjecture that our approximation for the wave 
functions does not alter order estimates of renormalons. 
In these approximations, the potential energy is given by 

d'^p r d'^k 



-^pot,/3o 



2M J (27r)4 J {27T 



xfkp)K,M xT^p + k) 



(0), 



Its Borel transform reads 



-^pot,/3oPJ 



/(n. A) = -2' 



m 



2u 



d^p 
(27r)3 J (27r) 



d^k 



(9) 



(10) 



2 + l)^ (^\^+ k\2 + [-(A;0A)2 + |fc|2 - ie 



+ 



.A;0 + l+p2/2+ \p + k\^/2-ie -k^ + l+p/2 + \p + k\'^ /2 - ie , 



(11) 



Note that in the integral of /(m. A), we have rescaled the 3-momenta p, k and the energy k^ 
by the Bohr scale pb = CFas{fJ')m/2 and the Coulomb binding energy Eb = (CFa5(/i))2m/4, 
respectively (j9° integration has already been performed). We have introduced a dimensionless 
parameter A = Eb/pb = Cfois/2, which characterizes the off-shellness of Q and Q in the 
boundstate. The Borel parameter u is defined with respect to /3oa5(/i)/(47r), so that the 
potential energy is given by 



E, 



(15) 
pot,/3o 



y , 

71=0 \ / 



-E. 



(IS) 



QyU pot,/3o 



\U\ 



(12) 



u=0 



If we take the limit A 
be evaluated easily as 



J{u) 



28+2«^^2 



in /(n. A) [Eq. (|TT])] before we perform the integration, it can 



d^k 



d^p 

(27r)3 J (277)4 



p 



+ l) {\p + k\^ + l 



2+2u 



X 



,A;0 + 1 +p 2/2 + |p + A;|2/2 - ie -k^ + 1 +p V2 + \p + k\'^/2 - ie, 



1 + 2m 



cos TTM 



(13) 
(14) 
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Rn X 



2 



/y»2 y 5 

^ 27rA 



21 



+ 



91 



8 7rA3 ^ 24 7rA 



X 



1287 



429 



256 7r AS 



1467 

128 7rA3 256 TT A 



/ 1215 
V 1024 TT 



X° X 



6149 



8195 



12155 I 6149 I 26169 j 

A^ 1024 ttAS 5120 TT A3 1024 tt A 



Table 1: The residues of the poles at m = in E^^''p^[u], normaUzed by _(£££^ljl}.-^ ^ il/[CFas{fJ^)m). 



If we replace /(m, A) by J{u) in Eq. (p!0D, it reduces to the Borel transform of the potential 
energy calculated from the static QCD potential in similar approximations.|3 In this case, IR 
renormalon poles are located at m = |, |, |, . . . on the positive real axis. 

Using Feynman-parameter integrals, we have reduced I{u, A) to a one-parameter integral 
form including hypergeometric functions: 

(15) 

= ^ h h my) - (-., 1; I; (1^)^) . (16) 



As can be verified easily, the integral in Eq. (15) is regular at Reu > and exhibits a simple 



pole as M — i> 0. Combining with the prefactor, we see that E'^ot]i3oW] ^as simple poles only at 
the positive integers m = (A^ = 1, 2, 3, ■ ■ ■) in the region Rew > 0. The residue Rn of the 
pole at M = is a {2N — l)-th-degree polynomial of A~^ and includes only odd powers of A~^. 
First few residues R^ are listed in Table |1|. 

The asymptotic behavior of the perturbative expansion of E^^^ at large orders is controlled 
by the poles at m = 1, 2, 3, ■ ■ -. The leading contribution comes from the pole at u = 1 and is 
given by 



4:71 



n r 



5" / Ri 



du"- \u 



«!. (17) 



Jm=0 n=0 \ ^" 



n=0 

We may estimate the uncertainty originating from the leading renormalon pole at u = 1 as| 



This can be obtained most easily by taking the expectation value of the Borel transform of the QCD 
potential in the large-/3o approximation with respect to the IS' Coulomb wave function (in position-space). 

oo 

^ Given the asymptotic series A = nl a" with a small parameter a > 0, we consider the finite sum 

n— 

up to n = n^, Aopt — ^^n!a", to be an approximate value of A, where n-i, « 1/a is the order n at which 

the term of the asymptotic series becomes smallest. Since the size of the term barely changes in the range 
n G {rii, — yjnl^ rii, + y^^T), we may consider that the estimate by the finite sum has a truncation error of order 
5 Aopt ~ ^/nl X nv,! a"* ~ e""*. (See e.g. [|lj.) 
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Similarly, the uncertainty generated by the renormalon pole at m = may be estimated as 



Ko%o - Aqcd X ^ y<{fo + fial + ha's + --- + fN-i cf-^) . (19) 



The polynomial of as^jj) follows from the form of the residue Rn'- fi denotes the coefficient of 
f in {af 

2a / A \ 2fe+l 

/Aqcd 
asm / \ a^m 



a'g in (a|^ Rn). We may rearrange the above uncertainties for all A^ and rewrite them as 

rp(is) A f Aqcd V" ( Aqcd \ .„„s 

<^^pot,A, ~ Aqcd X — — X ) , (20) 



where a, 6 > 0. Correspondence to Eq. (|T^ can be made via A^ = a + 6 + 1. Since in 
the on-shell approximation, the pole at m = a + ^ in J{u) induces an uncertainty of order 
Aqcd X (Aqcd/(«5^))^", we see that each renormalon uncertainty is suppressed by factors 
(AQCD/(a|"^))2''+i due to off-shell effects. 

Let us compare the above result with the order estimates of renormalons by power counting 
in the case a|m ^ Aqcd- We apply the following counting rules to Eq. (H): 

~ Aqcd, / ~ a|m, |p1 ~ a^m. (21) 

Also, we regard the kernel as K^^^ik) ~ as{\/ —k"^) /\k\'^ and count (only) this coupling as 
as{\/—k'^) ~ [For example, we count the integral measures as d^p ~ alm"^, d^k ~ ^qcti-i 

the Coulomb wave function as (pcisip) ~ (asm)'^^"^, etc.] We expand the propagators and 
4'c,is{p + k) in terms of k^ and k in the integrand. Due to the symmetries under ^ — 
and k ^ —k, odd powers of k^ or k vanish. Then the expansion generates power series of 
(Aqcd/(«|'^))^ and (AQCD/(tts'm))^ in the power counting. The order estimates of renor- 
malons agree with Eq. (^). 

In the opposite case where the off-shellness is small, a;|m ^ Aqcd, the power counting of 



|2 



^^QCD) 



Eq. (^ leads to a different order estimate. Since ~ -LL or a|m, and (fc°)^ — \k\ 

it follows that \k\ ~ Aqcd ^ 1^°! ~ or a|m (dominant contribution comes from an 

instantaneous gluon), and the leading renormalon is estimated to be of order Aqcd- Despite 
of this order estimate, the leading pole in Ep^^^^^[u] is located at m = 1 even for small off- 
shellness, < A ^ 1, and is not shifted to m = 1/2. To verify consistency, we now examine the 
uncertainty of the perturbative series of the potential energy when the off-shellness is small. 

First, by analyzing I{u, A) in the vicinity of A = 0, we find that it can be decomposed into 
the parts which are analytic and non- analytic at A = 0:| 

/(m,A) =/a(m,A) + /na(m,A). (22) 

The analytic part coincides with J{u) at A = 0: 

/a(m,A) = J(m) + 0(A). (23) 

^ It is well known that asymptotic expansions of loop integrals of Feynman diagrams consist of analytic and 
non-analytic parts in general. See e.g. |22]. 
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For u > 0, the non- analytic part behaves as 



/na(h, A) = A^-^" [ 4l{u) + 0{A) \ ( « > ), (24) 

where 

^0) . . _ 2r(| + 2..)r(i-^) 

J(ii) and A^~^"7]lf°2('^)) respectively, have the leading poles at = 1/2 on the positive real 
axis. These poles cancel with each other (independently of the value of A). By contrast, in 
the vicinity of m = 0, and if A <C 1, the contribution of the non-analytic term is suppressed 
and |/a| ^ I-^naI holds. 

Let us write the contribution of the non-analytic term I^)^{u) to Ep^^f^^[u] (up to a normal- 
ization factor) as 



X 



The asymptotic behavior of the expansion coefficient c„(A, x) for large n is determined by the 
pole and its residue at = 1/2 in the left-hand-side: Cn(A, x) ~ —2'"^'^x/7r. As A becomes 
smaller, c„(A,a;) approaches this asymptotic behavior more slowly. Stating more precisely, 
the asymptotic form is a good approximation if n ^ nx/A = log(a;/A). When we estimate the 
uncertainty of the perturbative expansion of the potential energy i^p^^ ^^, the size of the term in 
the vicinity of n = n^,, where it becomes minimal, matters. If ^ n^/^ (i.e. a|m ^ Aqcd), 
there are contributions from both /a and Jna around n = n*; the poles at n = 1/2 cancel and 
the estimate of I{u, A) by the renormalon pole at it = 1 is valid. Oppositely, if <^ Ux/a (i-e. 
Q;|m <^ Aqcd), the contribution from 7na is suppressed around n — n^, and the estimate by 
the u— 1/2 pole of Ia{u, A) J{u) is vahd. Thus, the behavior of the perturbative expansion 
at n ~ varies as if the position of the leading renormalon pole is shifted, depending on the 
relative magnitude of Ogrri and Aqcd- 

As we have seen, the leading renormalon contained in the potential energy is suppressed 
when the off-shellness is large. Then, an important question is whether the cancellation between 
the leading renormalons contained in the potential energy ii^pot and in the self-energies £'se still 
takes place in this case. Furthermore, if the cancellation occurs: Are the remaining sub-leading 
renormalons in the total energy also suppressed as compared to those in 2mpoie + Vqcd('")? We 

f 15) 

consider these questions. Using approximations consistent with those used for E'pQ^ jjj^, we may 
estimate IR renormalons in the self-energy contribution to the IS energy level as 



EsE,(3o -M-i -— \<PcMp) I M , ^0 — ~ p2 , . + ^ — - — -T- 



X 



2 ' ^ 2m ' 2 ^ 2m 

f+pO-m-£ - E(p) -Sm] [f - P° - m - £ - E(-p) - 5m 

2 



i?2 



^ 2 m + ('^^"'^^^^) ^^2Sm + + 0((E + Smf) , (27) 
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where 



S(p) 




X— — E(^^l (28) 

n=0 



represents the contribution of IR renormalons to the quark self-energy. The first two terms of 
Eq. (|27|) are the tree-graph contributions; 6m in the third term represents a scheme- dependent 
mass counter-term (which is free from IR renormalons for a short- distance mass); the fourth 
term represents the contribution from the self-energies of Q and Q [linear term in S] induced 

(IS) 

by Coulomb gluon in the large-/3o approximation. The Borel transform of reads 



2u 



with 

7 (27r)^7 (27r)4 [^-2 + 1^ ^ _(p^)2 + |^|2 _ 



1 



X 



+ 1+^2/2+ |p + A;|V2-ie -A;0 + l+pV2 + |p + A;|V2-ze, 

23+2" r(f + m) r(2 + 2 u) r(i - u) 



(30) 



3 7r3/2Ar(M)r(3 + u) 



X 



(1 - ' (1 + 0"'""' 2F1 (I + 1 - «; 3 + ^z; 1 - i,) . (31) 



It is easy to calculate the residue of at u = N {N = 1, 2, 3, ■ ■ ■) using the last line: 

the residue is a {2N — l)-th-degree polynomial of and includes only odd powers of 
(similarly to the residue of E^(^^i3q[u])- One may confirm that the order A^"^^ term of each 
residue is equal in magnitude and opposite in sign to the corresponding term of the residue of 
-^pot,/3oM at M = A^. Thus, the renormalons at order Aqcd x (^qcd/(«|'^))^^~^, which are 
contained in the potential energy and the self-energy contribution, respectively, get cancelled in 
the energy level M^'^^ = Eg^j^^^ + E^^^^^^^. This means that all the uncertainties corresponding 
to a = in Eq. (|20|) are cancelled. The leading renormalon contained in M^^ is of order 
Aqcd X (Aqcd/(«5"^))2 x (Aqcd/(«|"^)). Namely, renormalons contained in the energy level 
(when expressed in terms of a short- distance mass) are suppressed by a factor Aqcd/(«|'^) 
as compared to those in 2mpoie + Vqcd('^)- The orders of renormalons in M^J"^-* are consistent 
with estimates based on power counting of the momentum-space integrals (|lT]) and (|30|). 

It is necessary to investigate significance of the off-shell effects for the realistic heavy 
quarkonium states. Let us consider the IS' states of the (remnant of) toponium, bottomo- 
nium and charmonium, and set the quark masses m as 175 GeV, 4.73 GeV(= Mx(i5)/2), 
1.55 GeV(= Mj/^/2), respectively. Also, we choose the scales as 50 GeV, 2.49 GeV and 
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Order (PqC^s)"' term of E^ot'/So (LO in A-expansion) 


A 


n = 


n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


tt 


A 

A+NA 


-2.63 
-2.34 


-0.97 
-0.79 


-0.46 
-0.32 


-0.26 
-0.15 


-0.18 
-0.079 


-0.15 
-0.046 


0.087 


hi 


A 

A+NA 


-0.32 
-0.24 


-0.25 
-0.17 


-0.26 
-0.14 


-0.34 
-0.14 


-0.53 
-0.15 


-0.99 
-0.19 


0.18 


cc 


A 

A+NA 


-0.28 
-0.17 


-0.35 
-0.20 


-0.61 
-0.25 


-1.35 
-0.41 


-3.67 
-0.66 


-12.0 
-1.55 


0.30 



Table 2: The series expansion of the potential energy E^p^^^pg [Eq. (|T^)] when I{u, A) is replaced by J{u) [A] 
and by J{u) + A^-^'^I^liu) [A+NA]. We choose ^i = 50, 2.49, 1.07 GeV, m = 175,4.73, 1.55 GeV, = 5, 4, 3 
and ai^'\fi) — 0.130, 0.274, 0.448, respectively, for the toponium, bottomonium and charmonium 15 states. 



1.07 GeV, respectively.!]] Then the respective binding energies Eb = {Cpas{fi)yrn/4: are 
evaluated as 1.3 GeV, 0.16 GeV and 0.14 GeV. We may consider them as typical sizes of 
off-shellness and compare them with Aqcd = 0.2-0.3 GeV. It seems certain that the off-shell 
effects are important in the toponium energy level. As for the bottomonium and charmonium 
states, more detailed analyses may be needed to clarify the significance of the effects. Here, 

(IS) 

as a simple analysis, we examine the series expansion of the potential energy E^^^ '^^^ [Eq. (|T2p] 
when we replace /(m. A) by J{u) (on-shell approximation) and by J{u) + A^~^"/na(w)- The 
terms up to n = 5 are shown in Table |^. Within this approximation, we find that the off-shell 
effects cannot be neglected for all these states. One should note, however, that, since A are 
not particularly small for the bottomonium and charmonium states (see the same table), errors 
due to neglecting higher powers of A are estimated naively to be 20-30% for these states. 

For a systematic analysis of the realistic quarkonium energy levels, it is indispensable to 
calculate explicitly the first several terms of the perturbative series of E^^^^^^ and E^^j^^. The 
integral expression Eqs. (|l^) and (^) is not very suited for calculating the expansion coefficients 
of Ep^^'i^^[u]. We will report the calculation of these coefficients and their systematic analysis 
(together with the details of the calculations of the present paper) in a future publication. 

It has been known that the positions of renormalon poles are different between the observ- 
ables constructed from the on-shell and off-shell Green functions [Q], and that the orders of 
renormalons can be different in perturbative QCD and in effective theories [Q. In this paper 



we have shown, through an explicit calculation, how the orders of renormalons are "shifted" in 
the series expansion of the potential energy when the off-shellness of Q and Q is incorporated 
and varied. We find the mechanism fairly interesting. As far as we understand, all the previous 
perturbative analyses of heavy quarkonium states do not incorporate the off-shell suppression 
effects. This may originate from a failure to identify all terms of the asymptotic expansion in 



^ The scales for the bottomonium and charmonium IS* states are taken from g]; the scale for the toponium 
IS state is taken as (roughly) the scale where the perturbative expansion becomes most convergent from the 
analyses of p3|. 
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A of Eqs. (^) and (|2^) , in particular the non-analytic terms (cf. Ref. expand the 



integrand of Eq. (pAf ) in A before integration, some terms of the asymptotic expansion are lost. 
We find that cancellation of renormalons in the total energy -E'pot,/3o + -^se,/3o is realized at a 
deeper level in perturbative QCD than what has been known previously when the off-shellness 
is large. The remaining renormalon is of order Aqcd x {-^qcd/ {.C(s^)Y x {-^qcd / (ct's^)) ■ 

The dimension of this renormalon uncertainty is same as that of the non-perturbative cor- 
rection by the local gluon condensate of the Voloshin-Leutweyler type We note that a 
priori the two contributions are of qualitatively different nature. The relation between renor- 
malon uncertainties of perturbative computations and non-perturbative contributions in oper- 
ator product expansion is studied in detail for the i?-ratio As for the quarkonium energy 
level, Ref. Jl^ discusses the relation between the order Aqqj-, renormalon uncertainty in the 
static QCD potential and (ultraviolet sensitivity of) a non-local gluon condensate within the 
pNRQCD framework; it concludes that the former can be absorbed into the latter (the pertur- 
bative renormalon is cancelled against "the renormalon ambiguity of some non-perturbative 
effects"); it does not calculate the expansion of the matrix element of the gluon condensate in 
as, however. 
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